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Abstract.- In this paper we prove the existence of a large class of periodic 
solutions of the Vlasov-Poisson in one space dimension that decay exponentially 
as i — > oo. The exponential decay is well known for the linearized version of the 
Landau damping problem and it has been proved in ^4^ for a class of solutions of 
the Vlasov-Poisson system that behaves asymptotically as free streaming solu- 
tions and are sufficiently flat in the space of velocities. The results in this paper 
enlarge the class of possible asymptotic limits, replacing the flatness condition 
in [4| by a stability condition for the linearized problem. 
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1 Introduction. 

Landau damping is a remarkable property of collisionless plasmas. This effect 
was discovered in [9] and it consists in the exponential damping of charge oscil- 
lations in the plasma due to the combined effect of the electrical fields generated 
by the charges and the dispersion in the particle velocities. 

In more mathematical terms, this effect is usually studied using the Vlasov- 
Poisson system for negatively charged particles moving in a constant background 
of positive charges that makes the whole system electrically neutral. A linearized 
Vlasov-Poisson problem is then derived assuming that the inhomogeneities of the 
charge distribution are small. For rather general initial distributions of particles 
and initial velocities the corresponding charge disturbances of this linear system 
decay exponentially fast. 

The mathematical theory of the Vlasov-Poisson equation, including global 
existence results for general initial data in space dimension < 3 has been 
established in several papers (cf. [B], [T^, [TT]). 

Landau damping has been also extensively studied in the mathematical and 
physical literature (cf. [3], [5], [7], [H], [H], [IS])- The results in these papers hold 
only for the linearized version of the problem. The derivation of the exponential 
decay for the linearized Vlasov-Poisson system, in the cases where such a decay 
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takes place, relies heavily on the analytic properties of the initial data (cf. [5], 
[8], [H]). However, it has been proved in [13] that exponential decay cannot be 
expected, even for the linearized problem, for nonanalytic initial data. 

There are also physical arguments that explain why the energy of the charge 
waves tends to be transferred to the distribution of velocities of the particles 
(cf. [3], 0), yielding exponential decay of the charge disturbances. In these 
arguments the linearity of the problem plays an essential role, in particular 
it is possible to decompose an arbitrary disturbance in monochromatic waves 
using Fourier transform. However, it is not obvious how to generalize these 
physical arguments to the whole nonlinear problem, even for small perturbations 
of the linear case. Nevertheless, the numerical simulations in [TB] indicate that 
exponential decay of the charge disturbances might be expected for the whole 
nonlinear Vlasov-Poisson model and a large class of analytic initial data. 

A rigorous mathematical proof of the existence of solutions to the Vlasov- 
Poisson system in a circle for which the electrical field E decreases exponentially 
has been obtained in [1]. However, a general stability condition ensuring such 
exponential decay has not been obtained yet. 

We now precise the mathematical problem that we will be considered in this 
paper. Assuming that there is a background of charge uq to make the system 
electrically neutral, the mathematical problem that is usually considered in the 
study of Landau damping is the following one: 

ft + 'vfx + Efy ^ , -oo < w < cx) , 0<x<L , t>0 (1.1) 
f{0,v,t)^fiL,v,t) (1.2) 

/OO 

f (x^v^t) dv ~ uq , < X < L , t > (1.3) 
-oo 

^{0,t) = <f>iL,t) , 0,(O,t) = (/.,(L,t) , t>0 (1.4) 

E = (l)^ , 0<a;<L , t>0 (1.5) 

f {x,v,t) = fo{x,v) = fe{v) + go{x,v) , t = (1.6) 

where fe {v) satisfies: 

/>oo 

fe{v)dv = no (1.7) 



The physical basis on assuming the existence of the background of charge 
uq is that the system contains two types of particles with very different masses. 
The lighter ones that can move easily are described by means of the distribution 
f {x,v,t) . On the other hand, there are some heavier "ions" that cannot move 
so easily and are replaced by the charge background uq. 

Without loss of generality we can assume, using suitable units; 

710 = 1 , L^2n (1.8) 

Notice that (|1.7p implies that fe is a stationary solution of (|l.ip - (|1.5p . More- 
over, p.6p and (|1.7p imply: 

godxdv ^ (1.9) 
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In this paper we will obtain a large class of periodic solutions of the nonlin- 
ear Vlasov-Poisson system (|l.ll) - (|1.6p in a bounded domain for which the cor- 
responding charge density p and electric field E decay exponentially as t — > oo. 
The result is a perturbative one, in the sense that the derived solutions will be 
close to solutions of the linearized equation. The results are close in spirit to 
the ones in [J]. The main difference between the results in this paper and the 
results in [4] is that we replace the condition on the limit distribution in Lemma 
3.1 by a stability condition on the limit distribution that is equivalent to the 
stability conditions obtained in studies in the physical literature on the Landau 
damping problem. 

The main difficulty of the problem under consideration, even with the above 
mentioned smallness conditions, is that a naive linearization near an equilib- 
rium distribution fe (v) , say a maxwellian distribution, does not allow to show 
that the remaining nonlinear terms are small. This problem will be explained 
in detail in Section [21 The rationale behind this problem is that, although the 
electrical field and the charge density converge to zero, the distribution func- 
tion f {x,v,t) for the particles does not converge necessarily to the equilibrium 
distribution /g (v) . On physical grounds this might not be expected in the ab- 
sence of collisions. Actually, the dissipation of the energy contained in the field 
E {x, t) must result in the gain of kinetic energy of the particles of the system, 
as it was noticed in [3] (cf. also [5] for a clear explanation of this). However, in 
the absence of dissipative mechanisms there is no reason to expect for the long 
time distribution of particle velocities to approximate to the equilibrium distri- 
bution fe {v) or even to be spatially homogeneous. In fact, if the field E{x,t) 
vanishes fast enough as t —^ oo, the only restriction that we can expect for the 
long time asymptotics of / {x, v, t) is to behave like a "free streaming" function 
/oo {x - vt, v) . 

The previous discussion suggests not to linearize around the equilibrium dis- 
tribution, but near a free streaming function {x — vt, v) . It turns out that 
under suitable analyticity assumptions on fac and suitable smallness conditions 
ensuring that /oo is close to a stable equilibrium /g (v) , it is possible to ob- 
tain solutions of the Vlasov-Poisson system defined in t e (0, oo) , such that 
/ {x, V, t) ~ /oo {x ~ vt, v) as t ^ oo, and E {x,t) ^ as t ^ oo. 

It is interesting to notice that the linearization near the "free streaming" 
distribution /oo {x — vt, v) has some peculiar features. Indeed, the charge den- 
sity Poo {x, t) = J /oo {x — vt, v) dv — 1 is not identically zero for any t > 0, as 
it could be expected from the fact that / ^ /oo and ^ as t ^ oo. It turns 
out that Poo (^j ^) approaches to zero exponentially fast as t — > oo under suit- 
able analyticity assumptions on /oo that will be made precise later. The decay 
to zero of the density for the free streaming problem, due to the dispersion of 
the velocities, is a well known fact (cf. 2J). The argument in this paper indi- 
cates that this convergence to zero plays a crucial role in the Landau damping 
problem. However, the averaging property is not the only ground explaining 
Landau damping, since some additional stability conditions must be satisfied 
by the distribution fe. The conditions on fe that are required to obtain Landau 
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damping have been explained by means of intuitive physical arguments (cf. [5] , 
[8]). However, the role played by the averaging of the asymptotic distribution 
of particle velocities has not been considered, to our knowledge, in the physical 
study of this problem. 

The main result of this paper is the following: 

Theorem 1 // /e (v) satisfies a suitable stability condition for a linearized ver- 
sion of the Vlasov-Poisson problem there exist initial data /o (x, v) — fe (w) + 
go {x^ v) with go satisfying lll.9\) such that the corresponding solution of the sys- 
tem Jp - J775|) satisfies: 

\E{x,t)\ + \p{x,t)\<Ce-^* , 0<i<oo 

for some 7 > 0, C > 0. 

More precise conditions on fe, go will be given in Section [51 and Theorem 

The plan of the paper is the following. We recall in Section 2 some general 
properties of the classical linearized Vlasov-Poisson system. In particular, the 
difficulties that arise in trying to handle nonlinear terms with this linearization 
procedure will be explained. Section 3 describes formally the new linearization 
procedure. Section 4 formulates in detail a set of conditions for /e, goo that 
will play a crucial role in the rest of the paper. Section 5 analyses in detail 
mathematical properties of the linearized operator. Section 6 contains the proof 
of exponential decay for the nonlinear Vlasov-Poisson problem, in particular the 
proof of Theorem [1] 



2 The classical Landau linearized problem. 

In this Section we describe the linearized problem that was originally considered 
by Vlasov (cf. [T2]) and Landau (cf. 9J) and that has been studied also in the 
mathematical and physical literature with different levels of rigor [3], [S], [5], 
[nil [HI)- The goal is to explain shortly the kind of difficulties that arise in 
trying to derive exponential decay for the field E associated to the nonlinear 
Vlasov-Poisson problem taking as a starting point the linearized problem. More 
precisely, we will explain why a naive linearization argument does not allow to 
derive exponential decay for the field. A linearization procedure that allows 
to obtain such exponential decay estimates for the nonlinear problem will be 
introduced in Section [3] The procedure that we will use to solve the linearized 
problem is equivalent to the one used by Landau, although it is different, but 
closer to the approach that we will use in the forthcoming sections to study the 
nonlinear problem. 
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We write: 



f{x,v,t) = fe {v) +g{x,v,t) 



We assume that 5 is a small perturbation of the steady solution /e (v) . Neglect- 
ing formally quadratic terms in g we obtain the following problem that is the 
one usually considered in the study of the linearized Landau damping problem: 



9t + vgx + Efe,v = , 0<a;<L , t>0 , v eR 

<Pxx^ I gdv , E ^<j)^ , < x < L , t > 



(2.1) 

(2.2) 

(2.3) 
(2.4) 
(2.5) 



g [x, V, t) = goix - vt,v) - / E (x - v {t - s) , s) fe,v (v) ds 

Jo 

where we assume that go, E are extended periodically to x € M with period L. 

Plugging this identity into (|2.2p we obtain the following integro-differential 
equation for the electric field: 



g {x,v,0) = go {x,v) , 0<x<L , u G M 

(/.(O,t) = 0(L,t) , 0^(O,i) = 0,(L,t) , t>0 
g{Q,v,t) ^ g{L,v,t) , -y e R , t>0 

Then, using Duhamcl's principle to solve (|2.ip . (|2.3[) we obtain: 



Ex {x,t) = / g{x,v,t)dv 



go (x — vt, v) dv- 

E {x — V {t — s) , s) fe^v (v) dv ds 

It is easier to study this equation using Fourier variables: 

00 

E{x,t)^ fcn(<)e™^ 

n— — 00 
00 

go (x, -y) = Y gn (v) e"^ 

n— — 00 

where the periodicity of (j) requires bo (t) = 0. 
Then the linearized problem becomes: 



(2.6) 



inbn (t) = 
We write: 

Then: 



dvgn (v) e-™"* - / 6„ (s) 
Jo 



^(0 



e-''«Ve,. {v)dv 



ds 



(2.7) 



5 



nt poo 

inbn{t)+ / bn (s) K {n {t ^ s j) ds = Gn (t) = / dw5„ (w) e""^* (2.8) 

Jo J -co 

In order to solve this equation we introduce a family of fundamental solutions 
given by: 

inBnit;to)+ f BrAs;to)K{n(t-s))ds^S{t~'to) (2.9) 
Jo 

where to > 0. It then follows that: 

poo 

hn{t) = / Bn{t;s)Gn{s)ds (2.10) 



This family of convolution integral equations (|2.9p can be solved using Laplace 
transforms. Given a function tp (t) in {t > 0} we define: 



(z) = / if (t) e-^^dt 
Jo 

We recall that the Laplace transform of the convolution: 

{ip * -0) (i) = / ^{t ~ (s) ds 
Jo 

is given by: 

{(p*tp) (t) = (^(z)V'(z) 

On the other hand the Laplace transforms of the sequence of functions 
Kn (t) = K (nt) are given by: 



Kn (t) = / K (nt) e-'^'dt = -k(- 
Jo n \n 



Taking the Laplace transform of (j2.9p it then follows that: 
inlfn {z; to) + -K (-) Bn (z; to) = e""*" 



Then: 



Bn(z;to) = - r^TTY (2.11) 



Notice that the convolution structure of (|2.9p implies that i?„ (i; s) — i?„ [t — s; 0) . 
This can be seen also using (j2.1ip . Moreover, classical properties of analytic 
functions imply that the absence of zeros of Qn {z) = in + (;|) in the half- 
plane {Re (z) > —7} for some 7 > implies the exponential decay of Bn {t; 0) as 
t — > 00. The function Qn (z) is, up to linear changes of variables, the so-called 
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Landau function, has been studied extensively in several papers. Notice that an 
equivalent way of writing it, using (|2.7[) . is the following one: 



Jo J-oo J-oo Z + IV J_ 



fe jv) dv 

-co (z + iv)^ 



The following results can be derived from the results in the paper (cf. 
also [5] for related, although more formal results). 

Theorem 2 Suppose that the function fe (v) is analytic in the strip |Im(w)| < 
A, go {x,v) is analytic in |Im(u)| < A, |Im(a;)| < A and satisfy 

\fe{v)\<Y^-^ , a>l , \lva{v)\<A 

\goix,v)\ < ^ , a>l , \lm{v)\<A, \Im{x)\<A 

1 + \v\ 

for some A, B > 0. Assume also that the function Qn (z) does not have zeros in 
the half-plane Im (z) > for some value of n. Then the n-th Fourier coefficient 
of E {x,t) defined in V2. 6\) decreases exponentially as t ^ oo. 



Theorem 3 Suppose that f^ satisfy the assumptions in Theorem\^ Assume 
that the function Qn (z) has a zero in the half-plane Im (z) > for some n € Z. 
Then, there exists go{x,v) satisfying the assumptions in Theorem\^ such that: 

for some 7 > 0, C > 0. 



Remark 4 Notice that the assumption on the zeros of Qn {z) in Theorem\^is 
not satisfied for all the nonnegative initial distributions fe (v) . There are several 
examples of such functions yielding instabilities in Jl^ . For instance, an 
example of function satisfying the analyticity assumptions on Theorem and 
where the corresponding function Qn {z) has zeros in the half-plane Im (z) > 
forn^Q is (cf fl^): 

There exist distributions fe (f ) for which the corresponding functions Qn (z) 
yield stability. For instance, the maxwellian distribution fe {v) = ^ has been 
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extensively studied in the physical literature and it has been shown to satisfy the 
assumption in Theorem\^ (of. JSjj , On the other hand there is a proof in 

of the fact that all the distributions of fe (v) with only one maximum yield 
stability for the corresponding linearized problem. 

An example, not so relevant physically, but where the condition on the roots 
ofQn{z) in Theorem\^can be checked easily is: 



fe[v) 



In this case, the function Qn (z) can be computed using residues: 



Qn {z) = i 
the zeroes of Qn (z) are: 



1 

n H 



(f + .gn(Re(f))) J 



z = — n ± i 

n 



The following result can be obtained using the inversion formula for the 
Laplace transform, as well as the methods in [T2] . 

Theorem 5 Suppose that /e satisfies the assumptions in Theorem and the 
function Qn (z) does not have zeroes in Ini(z) > 0. Then the function _B„ {t;to) 
defined by means of i2.9\) satisfies: 

\Bn{t;to)\<-e-''-*-'°^ , t>to ,B„(i;to) = , t < to 
n 

This estimate for the fundamental solution i?„ (t; to) implies an exponential 
decay for E and for initial data go satisfying the assumptions in Theorem O 
Indeed, these assumptions on go imply: 

< C„ee-^l"l* 

where the constants C„ decrease faster than any power law if the function 170 is 
assumed to be C°° in x. 
Then, using (PTU]) : 

\bnit)\<Cnee-'' 
for some 6 > 0, whence the exponential decay of E follows. 

Let us now indicate the difficulty that arises in trying to use the results for 
the linearized problem in order to derive decay of the solutions for the whole 
nonlinear problem. Suppose that we keep the neglected quadratic terms in the 
linearized equation (|2.ip . Then: 

9t + vg^ + Efe,v + Eg,, = , 0<a;<L , t>0 , v (2.12) 



\Gnit)\ < 



dvgn [v) e 
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Using Duhamcrs formula it then follows that: 



{x, V, t) = go{x - vt,v) - / E{x - v{t~ s) , s) fe^v (v) ds 



— / E {x — V {t — s) , s) g {v) ds 
Jo 

Plugging this identity into (|2.2p as before it follows that: 



Ex {x, t) = g {x, V, t) dv — I go {x — vt, v) dv 



E (x - V {t - s) ,s) fe,v (w) dv 



ds- 



E [x — V {t ~ s) , s) g {x, V, s) dv 



ds 



Notice that, due to the presence of the last term containing the nonlinear 
part of the integration Jp [...] ds it might not be possible to ensure that this term 
converges to zero exponentially fast. Actually, this does not happen in general. 
Therefore, it might not be possible to treat this term as a small quadratic per- 
turbation of the above considered problem. This difficulty is not resolved by 
using more sophisticated linearization procedures, for instance approximating 
the characteristics of the whole nonlinear problem by the free streaming char- 
acteristics plus a corrective term and expanding by using the Taylor series in 
powers of the corrective term. The reason behind the failure of these approx- 
imation procedures is that the function g{x,v,t) does not converge to zero, 
although the electric field E does. As a consequence, all these linearization 
procedures contain terms that may not be expected to be small as i — > cx3 and 
some method for handling them must be found. 



3 Linearization near infinity: Formal computa- 
tion. 

In this Section we describe at a formal level a more convenient way of linearizing 
the Vlasov-Poisson system in order to derive estimates for the nonlinear terms. 
The key idea is to linearize around the expected asymptotics of the solutions as 
t — *■ oo. However, the function / (x, v, t) is expected to be oscillatory as t —^ oo, 
since in the absence of the field its dynamics would be given by the transport 
equation, and then we might expect: 

f {x,v,t) ^ F (x — vt,v) as f — > cx) 

for some function F. In order to obtain a function converging to a limit as 
t — !■ cxD it is convenient to introduce a new set of variables: 
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z — X — vt (3-1) 
f{x^vt,v,t)^f{x,v,t) (3.2) 

Then: 

ft ft 7 fx fz ; fv fv ifz 

Then (|l.ip becomes: 

ft {z, V, t) -tE{z + vt, t) f, (z, v,t) + E{z + vt, t) f, (z, v,t) = {) 
On the other hand we can rewrite (11.31) as: 



Ex{x,t)= I dvf{x,v,t) — l= / f {x — wt,w,t) dw — 1 
where we have used that ng = 1. Also: 

/>oo 

Ez {z, t) = f {z — wt, w, t) dw — 1 



From now on, we will drop the tilde from / in order to simplify the notation. 
Therefore we need to study the problem: 



ft (z, V, t) -tE{z + vt, t) /, (z, v,t) + E{z + vt, t) fy (z, v,t) = (3.3) 

/OO 
f {z -wt,w,t)dw -1 (3.4) 
-OO 

More precisely, we will construct solutions / {z,v,t) satisfying: 

f {z,V,t) ^ fe{v) + gooiz,v) = foo{z,v) 8.8 t ^ CO (3.5) 

where the functions fe (v) , goo [x, v) will be assumed to satisfy suitable analyt- 
icity assumptions that will be made precise later, and the function 1700 (x, v) is 
assumed to be of order of the small parameter e and: 

gao{z,v)dz = (3.6) 





The characteristic equations associated to (|3.3p satisfy: 

Z (t,t; z,v) = z , V {t,t; z,v) — V (3.7) 

and: 

dZ dV 

— {s,t;z,v)^-sE{Z + Vs,s) , — {s,t;z,v) = E (Z + Vs,s) (3.8) 

OS OS 
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We define the functions: 

Z {oo,t; z,v) = Zaa {t; z,v) , V {oo,t; z,v) ^ Voc (t; z,v) (3.9) 
Tlien, the solution of p.3p can be written as: 

/ {z, V, t) = /oo (^oo (<; z, v) , Voo (t; Z, v)) 
Using (|3.5p it follows that: 

/ (z, W, t) = /e (Ko (t; Z, V)) + goo {Zoo (t; Z, v) , Voo {t; z, v)) 

Therefore: 

/oo 
[/ {z - wt, W, t) - fe {w)] dw 
-OO 

whence: 

/oo 
[fe {Voo {t; Z - wt, W)) - fe {w)] dw (3.10) 
-oo 

/oo 
ffoo (^oo {t; z ~ wt, w) , Voo {t; z - wt, w)) dw 
-oo 

In order to linearize (|3.10p we argue as follows. We are trying to construct 
solutions where E decreases exponentially as t ^ oo. Moreover, the field can 
be expected to be of order e. In particular the functions Z, V can be expected 
to be nearly constant. We can then approximate the solution of the equations 
(|3.8p to the leading order as: 

/oo 
sE{z + vs,s)ds (3.11) 

/oo 
E{z + vs,s)ds (3.12) 

Linearizing the functions /e,5oo using Taylor's theorem, and taking into 
account (|3.1ip . (|3.12p we obtain: 



/•oo poo poo 

Ez{z,t)— I ds fe^v{w)E{z — wt + ws,s)dw+ / goo {z — wt,w) dw 

J t J —oo J —oo 

poo I'OO 

— I sds I °° {z — wt, w) E {z — wt + ws, s) dw+ 

Jt J-oo (^Z 

r°° r°° dgoo 

+ ds °° (z — wt, w) E (z — wt + ws, s) dw (3.13) 

Jt J-oo 



Equation (|3.13p is the linearized problem as < ^ oo. Our goal is to show that 
this problem provides a good approximation for the solutions of the original 
problem (frT|) - (fr6|) as t —> oo. 
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It is relevant to notice that (|3.13p contains terms that can be expected to be 
of order namely the last two ones. On the other hand, the first two terms 
on the right hand side of (|3.13p can be expected to be of order e. The reason 
for keeping the last two terms is that due to the presence of the term s in the 
third term of p.l3p it is not "a priori" clear if this term can be neglected for 
times t w i or larger. The last term in p.l3p has been kept to ensure that for 
the resulting problem J^^ {z,t)dz — 0. Indeed, assuming that this identity 

is satisfied and taking into account that J^^ E (z, t)dz = Q as well as p.6|) it 
follows, after integrating (|3.13p va z € [0, 2tt\ : 



2iT I'oo coo 

dz I sds 

Jt J-oo 

It: 



dgo 



{z — wt, w) E [z — wt + ws, s) dw+ 



dz / ds 



dgo 



dw 



(z — wt, w) E (z ~ wt + ws, s) dw (3.14) 



In order to check this inequality we notice that: 

Dgoc , . ^ dgoo , , ^ , dgoo I . ^ 
— — — (z — wt, w) — ~t— — [z — wt, w) H — - — [z — wt, w) 
Dw oz ow 

Therefore, the right hand side of (I3.14p can be written as: 



(3.15) 



/'27r j*oo 

ds dz dwE {z — wt + ws, s) 

t Jo J-oc 



[t — s) — — [z — wt, w) H — — — (z — wt, w) 
oz Dw 



and, integrating by parts we can transform (j3.16p into: 



(3.16) 



ds dz dwgoD {z — wt, w) ■ 



[(t — s) Ez {z — wt -\- ws, s) — (< — s) Ez (z — wt + ws, s)] = 



whence (|3.14p holds. 

To conclude this section, we remark that there is a way of writing (|3.13p 
where it becomes apparent that the last two terms are really small perturbations. 
Indeed, using (|3.15p we can rewrite (|3.13p as: 

/•oo fOO /"OO 

Ez[z,t)— I ds fe,v{w)E{z — wt + ws,s)dw+ / goo [z — wt,w) dw 

J t J —oo J —oo 

noo 



-ds 

t 



t -J —oo 

oo pOO 

ds 

t 



Dgo 



Dw 

s\ dg, 



[z — wt, w) 



E (z — wt + ws, s) dw+ 

(l ^ (z — wt, w) E (z — wt + ws, s) dw 

V t/ ow 



12 



and, after integrating by parts in the third equation: 



Ez{z,t)— j ds I fe^v (w) E {z — wt + ws, s) dw + I goo (z — wt,w) dw 

(3.17) 



+ / ~ ^ids I goo {z — wt, w) Ez {z — wt + ws, s) dw+ 



oo 



ds [ (l ) {z — wt, w) E (z — wt + ws, s) dw 

J_oo ^ t/ dw 

Suppose now that we study (|3.17p in a space of functions satisfying \E {z, t)\ + 
\Ez {z, t)\ < Me-''* for some suitable 7 > 0, M > 0. Then, the third and fourth 
terms on the right side of p.l7p can be estimated as: 

It then follows that, for a suitable choice of 7, the last two terms in (|3.17p 
can be expected to be small perturbative terms. This fact will be made rigorous 
in Section [SI 

Therefore, we expect to be able to approximate (|3.17p as: 

/•oo /"OO /"OO 

Ez{z,t)— / ds ff..v{w)E{z — wt + ws,s)dw+ / g^o (z — wt,w) dw 

J t J —OQ J —OQ 

(3.18) 

Equation (|3.18[) is basically equivalent to the linearized Landau damping 
problem studied in Section [21 except for the fact that the linearization has been 
made at t = 00. We will study first some properties of the linearized problem 
(|3.18p using Fourier analysis. We will then study the whole nonlinear problem 
(|l.ip - (|1.6p using a suitable functional framework. 

4 Analyticity assumptions for /g and 5^00 and some 
consequences . 

We make precise now the assumptions made on the functions /e, g^o- 
Assumptions (A): 

• The function /e {v) is analytic in the strip |Im {v)\ < A and it satisfies in 
that set 

|/eWI<— %r , a>l , (4.1) 
1 + 

• The function g^o ix,v) is analytic in the sets |Im(a;)| < A, |Im(w)| < A 
and it satisfies in this set: 

\gooix,v)\ < ^ , a>l , a 7^ 2 (4.2) 
1 + \v\ 
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• The function goo {x, v) is periodic in the x variable with period 2tt and it 
satisfies: 

goo {x, v)dx = (4.3) 

The requirement a 7^ 2 seems at a first glance a bit artificial. This assump- 
tion has been made in order to avoid the onset of logarithmic terms that would 
introduce nonessential technical difhculties. Notice that we can always assume 
that a ^ 2 reducing the value of a a bit if needed. 

There is a consequence of (|4.2p that will be used repeatedly in the following. 



Lemma 6 Suppose that ^■2\ ), J^. j| ) are satisfied. Then the function 

/oo 
9oo (2 — wt, w) dw (4-4) 
-00 

satisfies: 

\H {z,t)\ <Cee-^* , z e [0, 27r] (4.5) 

where 7 > can be chosen arbitrarily close to A, j < A and C > depends 
only on a, 7, A. 

Proof. Due to the analyticity properties of the function goo we can rewrite 
H (2, t) for |Im (2)1 < A + At using contour deformation as: 



H {z,t) — I goo {z — wt,w) dw (4-6) 
Jr(z,t) 

where F (2, i) = {w e C : Im (w) — a (2, i)} with |Im (2) — ta (2, t)\ < A, 
\a{z,t)\ < A. It is readily seen that for any 2 satisfying |Im(2)| < A + At it 
is possible to choose such a (2, t) , although in a nonunique way. For instance, 
given 2 in |Im (2) — rt\ — 13 < A with |r| < A, |/3| < A we can choose a (2, t) = r. 
It then follows from the representation formula (|4.6p that H {z,t) is analytic in 
|Im(2)| < A + At. Moreover, using l|4.2p it follows that 

\H{z,t)\<Ce (4.7) 

in |Im(2)| < A + At. 

The periodicity of goo in 2 implies that H (z,t) is periodic in 2 for any t > 
and ()4.3p implies 

H{z,t) dz = 



We can then write H {z,t) using the following Fourier series: 



H{z,t) = Y.^ e-™«il(e,t)d^ (4.8) 
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Using the analyticity properties of (•, t) we can rewrite the Fourier coeffi- 
cients in (14.81) as: 



Jo J[0,27T]+iL 

where L e {—A — At, A + At) . Choosing sign (L) = —sign (n) and choosing 
L = J + jt with J < A arbitrarily close to A it then follows from (|4.7p that: 



r e-^<H{^,t)d^ 
Jo 



and plugging this estimate into (|4.8[) we obtain (14.51 



5 Construction of the fundamental solution of 
the linearized problem at t = oo. 



Equation ()3.18p as well as the form of the nonlinear terms that have been ne- 
glected in the derivation of it suggest to study the following problem: 

/>00 /"OO 

Ez{z,t)= / ds fe,v {w) E [z ~ wt + ws,s) dw + h{z,t) (5.1) 

Jt J -oo 

where h (z, t) is a bounded function decreasing sufficiently fast as t — s- oo and 
satisfying 



h (z, t) dz — 



(5.2) 



In order to study this problem we will construct a fundamental solution 
associated with it. More precisely we will derive, using Fourier analysis, an 
explicit formula for a solution of: 



/oo fOO 
ds 



fe,v (w) G {z — wt + ws, s) dw + 



5{z) 



1 

2^ 



satisfying: 



G'(z,t) = , t > , z e 
G{z,t) = G(z + 27r,t) 



(5.3) 



(5.4) 
(5.5) 



Using (|5.2p as well as the invariance of the homogeneous part of (|5.3p under 
spatial and time translations we obtain that a solution E {z,t) of (|5.ip can be 
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written as: 



ds i d^G{z-^,t- s)h{^,s) (5.6) 

-oo Jo 



/•oo /'27r 

/ ds d(,G{z-^,t-s)h{^,s) 
Jt Jo 



since, due to 



h{z,t) = / ds d^ 



2n r 1 



There is a function that plays a crucial role in the whole theory of Landau 
damping and that appears in slightly different forms in different papers devoted 
to this subject. This function, that is usually referred as the Landau function 
takes the following form in our setting: 



/ '^'"^ ' dw-n' (5.7) 
Jm w-rj 

This function is defined in e C\M , n — ±1, ±2, .... If /e satisfies (|4.1I) . the 
function (f> can be extended analytically to the domain {Im (ry) > —A} for any 
n = ±1, ±2, .... It is worth mentioning that the function <I> (77; n) is discontinuous 
for 7] € R. Indeed, due to the Plemej-Sokolski formula (cf. yy) : 

$ iVo + iO) - $ (770 - iO) = 27rz/e,„ (vq) 

Therefore the analytic continuation of the function $ to the domain {Im (77) > 
is not given by the integral formula (j5.7p . This is a well known property of the 
Landau function (cf. [8J). 

The main result of this Section is summarized in the following Theorem: 



Theorem 7 Suppose that satisfies (^TTP- < S < A. Define two func- 
tions ip^ (77) by means of: 

(^) ^ / l^^dw, (5.8) 

Jc±s "1^^ 

where C±s ^R±iS. 

We define a function Q (z, t) by means of: 

(5^(z,i) = ^ > / -. , ^ ^ drj (5.9) 

(27r)' ^ " ■''-00 fi^^fi^M^"! 

V ^ n=-oo , n^^O |„|„ y 

Q{z,t) = J^ Q,{^,t)d^ + ^ ^Q.{i,t)di (5.10) 
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as well as a tempered distribution G by means of: 



{z,t)^Q, iz,t) + 6it) 



J2 S{z + 2T:e) 



2tt 



G{z,t) 



£g._ {^,t)d^ + l-J^\G. {U)d^ 



(5.11) 
(5.12) 

Suppose that the Landau function defined in ( [5. 7| ) does not have zeroes in 
the half-plane {Im (77) > 0} . Then, the series defining the function in 115. 9\) 
converges for any z e K, i G M. The distributions G, Q are supported in the set 
t < and they are periodic in z with period lix. We have Q G (Mx (—00, 0)) . 
Moreover, the following estimates hold: 



\Q{z,t)\ + \Q, {z,t)\ < 



C 



(1 



< t < 1 



-TT < Z < TT 



|Q(z,t)| + |Q,(z,t)| < Ce-"^!*! , t>l , zG 



(5.13) 
(5.14) 



where a is as in |^.j[ ) and the constants C, a depend on B, a. 
The distribution G solves HS.SD-^STT 



In order to prove Theorem [7] we begin by deriving the analyticity properties 
of the functions i/'i as well as some basic estimates for them. 

Lemma 8 Suppose that f^ satisfies J^.Jp and let 6 be as in Theorem Then 
the functions ip^ (77) defined in fi5.8\] are analytic in {Im (77) > —5} and they 
satisfy the following estimates: 



\m 



1 

Gk 



, Im (r/) > 



|/3 



Im(77)>--, fc = l,2,. 



(5.15) 



(5.16) 



1 + 1771' 

where 1 < f3 — min {a, 2} — Eq with Sq > that might be chosen arbitrarily 
small, G depends only on a, B, 6, Eq and Gk depends on a, B, S, Eq-, k. 

Proof. The analyticity of t/^j. (7;) in {Im (rj) > —5} is just a consequence 
of the analyticity of the functions for each w G G±s in the half-plane 

Im (77) > —(5 as well as the assumption (|4.ip . 

To derive ((5T5)) we use the inequality: 



Then: 



77 ± w 



fe,v (w) dw 



< 



< 



B 



m Jc±s l?7± w| 1 



\dw\ 



(5.17) 
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Using the fact that J(j^^ fe,v (w) dw = and spUtting the integral on the 

right side of (|5.17p in the regions where |w| < and jwl > respectively, we 
obtain: 



B 



w 



B r 1 



\v\ Jc±sn{\w\>^^^} \r]±w\ |w;r 



- \dw\ + 
- \dw\ 



(5.18) 



We can estimate the first integral on the right side of (|5.18[) by a constant 
if a > 2 and as C \ri\ if a < 2. (Notice that a 2, and therefore the 
logarithmic case does not occur). On the other hand we can estimate the second 
integral on the right hand side of (j5.18p introducing the rescaling w = C- 
Then: 

^ MCI 



^n{KI>U 



1^ 



icr 



The integral in this formula is convergent, since a > 1. The main contribu- 



tion to this integral for large I77I is due to the region where 



is small. 



The closest distance between C and =FrT is of order -r-j. Therefore: 

^ ^\v\ \m 



n{lCI>H 



icr 



MCI <c,-[|iog(H)| + i] 



whence, choosing eo > arbitrarily small, (|5.15p follows. 

Estimate (|5.16p can be proved combining (|5.15p as well as the fact that 
the functions {rj) are analytic in {Im {rj) > —6} and the classical Cauchy's 
inequalities for analytic functions (cf. [T]). ■ 



The functions ipj. (rj) are closely related to the Landau function $ (cf. (|5.7p ). 
We reformulate some properties of the functions {rj) in terms of properties 
of the function $ that have been studied often in the literature (cf. for instance 

my 



Lemma 9 Suppose that /g satisfies Suppose that the Landau function 

4> (77; n) defined in ( [<5.7p does not have zeroes in the region {Im (77) > 0} for any 
n — ±1, ±2, .... Then, there exist i/q > and 6 > such that 



lgn{n) 



iv) 



\n\ n 



> 



(5.19) 



for Im (ry) > — i^o- 
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Proof. Using (|5.8p we have: 

1 -Y = 1 T / — ■ — ■ T-— '^w^ = 2^ {-Sign {n) rj; n) 

\n\n n'^ J(j^^ w + sign[n)rj 

(5.20) 

Notice that (|5.7p as weh as the fact that /e {w) takes real values for w e R 
implies: 

$(7);n) = (<i> (77; n)) (5.21) 

Therefore, if $ (r/; n) does not have zeroes for {Ini [rj) > 0} , it does not have 
zeroes for {Ini(77) < 0} either. Moreover, due to (|5.2ip we can restrict our 
attention to the case sign (n) = — 1. Due to (|4.ip we can extend analytically the 
function $ (ji; n) to the domain {Im (77) > —A} by means of the formula: 

9[ri;n)= / — ■ aw — n 

It then follows, due to (|4.ip that for \ri\ > p with p independent of n: 



— $ {—sign {n) 77; 7T,) 



1 

> - 

- 2 



On the other hand, since <i> (77; 77,) does not have zeroes for r\ G {Im (77) > 0}, 
\'r]\ < P it follows by continuity that the analytic extension of ^ycj) (7/; ti) to 
{Im (77) > —A\ does not have zeroes in the region —vq < Im (77) < 0, I77I < p 
for some vq > sufficiently small. Notice that vq can be chosen uniformly in 77,. 
Using (|5.20p the result follows. ■ 

Remark 10 The absence of zeros of the Landau function $ (77; n) in the half 
plane {Im (77) > 0} for any n = ±1, ±2, ... is precisely the condition required for 
the stability of the solutions of the linearized Landau problem studied in Section 
Therefore, the functions fe {v) yielding stability for the problem considered 
there yield also stability for the linearized problem considered in this Section. In 
particular, the examples of stability and instability in Remark^ are also valid 
for the linearized problem considered in this Section. 



We now derive some estimates on the Fourier coefficients in the series 
that will ensure the convergence of the series. 

From now on we will write by shortness ± instead of sign {n) . 



Lemma 11 Suppose that /g satisfies J[ ) and let S be as in Theorem^ Sup- 
pose that the Landau function $ (77; 77,) defined in ( [5. 7| ) does not have zeroes in 
the region {Im (7/) > 0} for any n = ±1, ±2, .... Then: 



e*l"l''*7/;± (77) 



drj 



< 



C-y 1 



teR, 7>0,7i^0 



(5.22) 
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where C-y depends on a, B, S, 7. Moreover, we have: 



and 



00 1 



(-^ \n\n 

where C depends on a, B, S, 7 



dr] ^0 for t > 



drj 



(5.23) 



(5.24) 



Proof. Notice that, for any t — 1, 2, ... we have: 

1 d' 



J\n\r)t _ 



{i\n\tf 977^ 
Then, integrating by parts we obtain: 



e*l"l''V± (r/) , (-1)' 
-drj ~ 



00 1 



(^|n|t)^ 



V'± (jy) 



d?7 (5.25) 



Using (|5.19p in Lemma [51 it then foUows from (|5.16p that: 



'0± (ry) 



1 - 



|n|n 



< 



1 + 



(5.26) 



with P as in Lemma 13 Combining (|5.25p . (|5.19p . (|5.26p we obtain 



g»|n|r,t^^ (77) 



00 1 - 



I n I n 



< 



{\n\tf 



teM, 717^0, ^ = 0,1,2, 



(5.27) 



Notice that for ^ = (jOT)) follows immediately from (jSTTS)) . Estimate ((O^ 
follows for 76 (£ — 1, ^) by interpolation. 

The identity (|5.23p can be obtained using the fact that in the half plane 



Im {rf) > where the fmiction 



is analytic and it is bounded as 



Ct 



we have also the estimate |e'l"l''*| f,-\n\iva(n)t < ^ r^j^^^^ ([5:231) follows by 
deforming the contour of integration (— cx),oo) to (—00,00) + iR with i? > 
and taking the limit i? — > 00. 

Finally, we prove (|5.24p using the identity: 



-00 1 



\n\n 



drj = 



V'± {ri) dr] - 



1 



n n 



(^± (^))' 



00 1 



\n\n 



drj (5.28) 
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The first term on the right hand side of (|5.28p can be computed using 
4'+ {i]) drj = lim / / iWi;J — ^-dwdrj ~ lim ( dwf^ „ {w) ( — ^ — 

oo J-R Jc±s '7 ± W i?.^oo Jci, ' J-R V±W 

Since /_"^ — log (^3^7;;^ is bounded for large R and w G C±s we can 
apply Lebesgue dominated convergence Theorem to obtain: 



V'i iv) dr] = dwfe^y (w) lim 



C±s 



log 



R±w 



-Txi \ dwfe,v (w) = 

(5.29) 



-R±w 

Using (EUl), (jgrg)) and ({05)1 in ([g?^ we obtain ((Oi)) . ■ 

We can now prove the convergence of the series defining [z, t) in (|5.9p . 

Lemma 12 Suppose that the assumptions of Theorem^ are satisfied. The se- 
ries on the right hand side of i5.9\} is convergent for any i g M. The function 
Qz defined by means of i5.9\) is identically zero for i > 0, Qz{-,t) G C°° (R) for 
any t < and Qz{-,t) G (R) for t = Q. Qz{-,t) defined by means of \5.9^ is 
periodic with period 2tt for any t £ M and it satisfies J^^ Qz {z,t) dz = 0. 

Proof. This Lemma is just a consequence of (|5.22p - ()5.24p . ■ 

Notice that (|5.22p - (|5.24p are not strong enough to derive uniform estimates 
for the function Qz {z, t) or its derivatives if t is close to zero. This is made in 
the following Lemma, where the self-similar structure of Qz {z,t) is derived. 

Lemma 13 Suppose that the assumptions of Theorem^ are satisfied. Then 
Qz{z,t) = f,(^^+R{z,t) , zeM,|z|<^, -l<i<0 
where: 

\R{z,t)\ + \Rz(z,t)\<C , zeM, ieM 
Proof. We can write 



°° e 



-00 1 



: -^^dr^= / e'l"l''V± (?7) rf?? + (i) (5.30) 



where: 



1 e'N"* {^^ 

[t) = -T-r- / —7 ^ /„-)X dr] 5.31) 

I ?1 1 71 



Using (|5T5l) . (IST^ . ([Qi]) we have: 
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\^n(t)\<^ (5.32) 

The term e*l"l''*^/'± {v) dr] can be computed explicitly. To this end we 
derive a more convenient formula for this integral. Using (|5.8p we obtain: 

lim / e*l"l''*7/;± (?7) rf?? 



R 



i\n\r]t f fe.v {w) 



lim / e'l"!"' / ' dwdTj 
R^°°J-R Jc±s 

r t-R gi\n\rit 

hm / dwfe,v {w) / — ; dry 

^-'O" Jc±, J~R V±W 

/R gi|n|i7t 
— ; dv 
-rV±w 

We just need to compute the integral for t < 0. This computation can be 
made by using residues. More precisely, we can replace the contour of integra- 
tion [-R,R] by [-R,R] U {z : z = Re'", 9 € [-7r,0]} . The contribution to the 
integral of the half-circle disappears as R ~* oo, Indeed, the integrand can be 
estimated as ^ if |Im {ri)\ < A is of order one, and as — - for Im (rj) < —A. 
Using the fact that the length of the half circle is bounded by CR, it follows that 
the contribution of the integral to the half-circle disappears taking the limits 
R oo, A — + oo. Thus: 

r-R pi\ri\r)t _ ^ ^ 



fit i\n\rit 

lim / df] = -27r?e=F'l"l"'* = -27rie" 



whence, using also the analyticity properties of fe- 



l"l''*?/'± iv) dv = -27ri / dwfe,v (w) e""""* (5.33) 



Using (|5.30p . (|5.33p we can rewrite (|5.9p as: 

Q,{z,t) =LJoiz,t)+uji{z,t) , t<0 (5.34) 



{27T? ^ " 
Notice that due to (|5.32p the function uii {z,t) satisfies: 

1^1 (z,t)| + 1^1,2 (z,i)| < C , zeR, teR (5.35) 
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On the other hand the function luq (z, t) can be expUcitly computed. Indeed, 
let us define the distribution: 



Tiz,t) = - 



E 



^in{z—wt) 



n— — oo , n^O 

We have, in the sense of distributions: 



dT , ,1 ^ 



^inz ^—inwt 



71— — oo , n^O 

whence, using the formula: 



oo ^ 

6{z-wt + 2'kI) - — 



T{z,t;w)= f T,{^,t;w)d^+^ f ^T, , t; w) 
Jo ^'^ Jo 

we obtain: 

oo / z 1 

T (z, t;w)= ^ [x iz-wt + 2TTi) - x {-wt + 2^^)] + - 

^ oo 

+ 7;- (wt - 27r^) [X (27r + 1) - wt) - X (27r£ - wt)] 



where x(') denotes the characteristic function whose support is (0,oo) 
Writing 

/oo 
dwfe,v (w) T {z, t; w) 
-oo 

it then follows that: 



Wo I 



,{z,t)= V / * /e,. Hdw+— V / * /e,„ (W) (wt - 2^^) 

whence, after some integration: 

oo 

^o{z,t)= ^ 



dw 



l=-oo 

z ^ „ /27r(^+l)\ t °° 



Using (|5.34p we obtain: 



- y 

27r ^ 

— oo t 



/e (w) dw 



Q,(z,t) = /e(|) +i?(z,t) 
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Riz,t) = -f,{0)+ 



z 

2^ 



fe 

271 (e + i) 



z + 27r£ 
t 



fe 



2t:£ 



t 



t 

2^ 



E 



Using dm]) and (jOS]) it follows that: 

\R{z,t)\ + \R,{z,t)\<C , ze 



Z < TT 



fe (w) dw + UJi {z,t) 



-l<t<0 



Thus the result follows. ■ 

We can now derive estimates for {z, t) for t large. 

Lemma 14 Suppose that the assumptions of Theorem hold. Then: 

\Q{z,t)\ + \Q, {z,t)\ < Ce'^' 

for < < — 1, with a > and C > 0. 

Proof. Combining (j5.15[) and (|5.19p we can write using contour deformation: 



drj 



ai\n\rit 



■0± (r)) 



-»7 I 1 



drj 



thus we have 



\n\n 



drj 



< Ce^l"l* , t < 



for some 7 > 0. It then follows from (|5.9|) that Qz {z,t) is periodic in z with 
period 27r, as well as analytic and bounded in |Im(z)| < -^t t < ^1- Arguing 
as in the proof of Lemma [S] we obtain that Qz {z, t) satisfies: 

\Qz{z,t)\<Ce''' ,t<-l 

for some a > 0. A similar estimate for Q (z, t) then follows from the formula 

Q(z,t)= r Q,{^^t)d^+^ f\Qz{Lt)d^ 
Jo Jo 

and the result follows. ■ 

To conclude the proof of Theorem [7] it only remains to show that the distri- 
bution G defined in (|5.1ip . (|5.12p solves (|5.3p - (|5.5p in the sense of distributions. 
To this end we will reformulate (|5.3p in term of the Fourier transform of G and 
we will verify that the solution of the resulting equation is the one given by the 
Fourier transform of G that can be computed using ()5.9p - (l5.12p . 



24 



Lemma 15 Suppose that the assumptions of Theorem^are satisfied. Then the 
distribution G{z,t) defined by means of HS.ll]) . i5.1S\) solves i f5. j^j - lJTS)) in the 
sense of distributions. 



Proof. We expand G {z, t) using Fourier series: 

n—oo 

G{z,t)^ ^ .9n(i)e» 



(5.36) 



Plugging (|5.36p into (|5.3p we obtain, after some computations that the func- 
tions gn {t) satisfy: 



ingn (t) 



/CO 1 
/e,„He-"™(*-^)d«; + — , n^Q (5.37) 
-no ^TT 



In order to solve these equations we compute the Fourier transform of the 
functions gn (t) that we define by means of: 

1 



?n (0) = 



We define also the functions: 



gn (t) e-^'^dt 



/OO 
V± it)' 
-OQ 



itd 



dt 



(5.38) 



where: 



f± it) 



fe,v{w)e^''"*dw , t<0 



ip^{t)^0 , < > 

The functions ip^. decrease exponentially as t — > — oo. This can be seen using 
(|4.ip and contour deformation to derive the following representation formula for 
these functions: 



^± it) = 



fe,v iw) e-'^^'dw , t<0 



(5.39) 



C±s 



where C±s ^ R ± iS, < S < A. 

Applying the Fourier transform to (|5.37p and using the fact that a convolu- 
tion is transformed to a product by this transformation we obtain: 



ingn iO) 



1 



7$^ 



9n ie) 



1 1 

27r 



(5.40) 



On the other hand, using (j5.39|) we can rewrite the functions <i>± in the 
following form: 



$± {0) 



-ite 



fe,v iw) e^^'^dw 



C±5 



dt 



fe,v iw) 



fe.v iw) 

c±. d±w 



dw 
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Then ([g^ implies: 



5n {0) 



1 



1 



(2^)' * n 



1 

|n| 



1 



or, using the inversion formula for the Fourier transform: 
1 



9n (t) = 



{2ti)iJ- 



1 

R 



-±W 



It is readily seen that this formula is the same that the one for the Fourier 
coefficients of the function G in (|5.12|) . Notice that the choice (|5.12p implies 
9o (t) = 0. This concludes the proof of the result. ■ 



End of the Proof of Theorem [Tl It is just a consequence of Lemmas 

laisi ■ 

Using Theorem [7] we can prove the following result that will play a crucial 
role in the analysis of the nonlinear problem: 

Proposition 16 Suppose that the assumptions of Theorem hold. Let h G 
C (Rx R+) be a function satisfying 15. h{z + 27r, t) = h (z, t) for z €R, t> 
0, as well as the estimate: 

\h{z,t)\ < Se"^* , z eR , t>0 

with < 7 < a. 

Then there exists a function E{z,t) G C (M x R+) solving h5.1\l , satisfying 
J^'^ E {z, t) dz^O, E{z + 27r, t) ^ E (z, t) and: 

\E{z,t)\ + \E,{z,t)\ <CBe-^' 

for some C > which depends only on A, B, a. 

Proof. The desired solution E {z,t) can be obtained by means of the formula 
(|5.6p . Notice that we just need the values of h {z, t) in t > to obtain E {z, t) 
for t > 0. Moreover, due to the linearity of the estimate we can assume B = 1. 
It is readily seen that J^^ E {z, t) dz = 0. Differentiating (|5.6p with respect to z 
we obtain: 

/oo /'27r 
ds / d^G,{z-^,t- s)h{^,s) 
-oo Jo 

Using (|5.1ip it then follows that: 

/oo /'27T 
ds / d^Q^{z-^,t- s)h{^,s) 
-oo Jo 



26 



Using (EUni), ([Oil it then follows that: 

\E,{z,t)\<Ce-'^' 
Finally we use the representation formula: 

Jo 27r Jq 

to obtain: 

\E{z,t)\ < Ce-^' 

whence the result follows. ■ 

6 Nonlinear problem: On the existence of expo- 
nentially decaying solutions. 

In this Section we will prove the existence of a large class of solutions of the 
nonlinear problem p.ip - (|1.6p that decrease exponentially fast as t — > cx). 

Theorem 17 Suppose that /e, goo satisfy Assumptions (A). Assume that the 
Landau function defined in |5. ?p does not have zeroes in the half-plane {Im (77) > 0} . 
Then there exist Eq = Sq (a, A, B) > 0, L = L (a, A, B) such that for any f^, g^^ 
satisfying Assumptions (A ) with £ < £0 there exists a solution f (z, v, t) of i3.3\) - 
\3.4\ l defined for z G M, u G K and < t < 00 and satisfying i3.5\) as well as: 

\E{x,t)\ < Cee-^^ for < t < 00 

Using the change of variables (|3.ip . (|3.2|) we can prove the following result: 

Corollary 18 There exists a function fo e (M x M) , with fg {x + 27r, v) = 
fo {x,v) such that the corresponding solutions of the system !il.l\} - !!77^) are de- 
fined for < t < 00 and they satisfy: 

\E {x, t) I < Cee-^^ for Q<t<oo 

Remark 19 In all the following C will denote a numerical constant depending 
only on a, B, A that might change from line to line. 

Proof of Theorem 1171 The strategy is to solve (I3.10p by means of a fixed 
point argument in a space of functions satisfying: 

\\E\\^snp{e'^'[\E[z,t)\ + \E,(zM] <^ 
t>o 

Suppose that < 00. Then, integrating the characteristic equations p.8p 
it follows that: 

\Voo{t;z,w)-w\<C\\E\\e-'^* , < i < 00 
\Zooit;z,w)'-z\<C\\E\\{t + l)e-''* , 0<t<oo 
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Using Taylor's expansion in (|3.10p we obtain: 



{z, t) = fe.v (w) (Voo {t; z - wt, w) ~w)dw + / goo {z - wt, w) dw+ 

+ [ {z — wt, w) {Zoo (t; z — wt, w) — (z ~ wt)) dw+ 

— — (z ~ wt, w) [Voo (t; z — wt, w) — w) dw + R [z, t) 
„o ow 

where: 

\R{z,t)\ <C||i;||'(t+l)'e-2^* 

In order to derive better approximations to the differences Zoo [t] z — wt, w) — 
[z — wt) , Voo (i; z — wt, w) — w we need to derive estimates for the solutions of 
the characteristic equations: 

dZ 

— = -sE{Z + Vs,s) (6.1) 

OS 

dV 

— = E{Z + Vs,s) (6.2) 
OS 

Z {t,t;z,v) ^ z , V {t,t;z,v) = V (6.3) 

with s > t. 

Using Taylor's expansion we can write: 

dZ ^ , . 

— = —sE [z + vs, s) 
ds 

+ o(^s \\E\\ e-^^ \Z-z\ + s' \\E\\ e""^^ \V - v\ + \\E\\ .e'^^ exp ("^e^* 

—— = E [z + vs,s) 
ds 

+ O (^\\E\\ e-T^ \Z-z\ + s \\E\\ e-^' \V - v\ + \\E\\ e'i' exp ("^e^* 

for some a > 0. The first two terms of the remainder arise in making a Taylor 
expansion for the range of values where sj^ — w| < 1. The last term is the 
contribution from the region where s |F — w| > 1. Since |y — w| < CM \\E\\ e~'^* 
this requires s huge, and due to the exponential decay of E the estimate follows. 
Taking into account the estimates for the differences \Z — z\ , \V — v\ that are 
similar to the estimates obtained for Zoo, Voo it then follows: 

Zoo{t;z,v) - z = - sE{z + vs,s)ds + 0(^{t+l)^\\Efe-'^'^^'^ (6.4) 

Vooit;z,v) - V = E{z + vs,s)ds + (^it+l)\\Ef e-^^*'^ (6.5) 



28 



We can then write the following equation for E (z,t) : 



OO /'OO 

ds 



sds 
ds 



fe,v (w) E (z + ws — wt, s) dw + I goo {z — wt, w) dw 
dgo 



— OO 

OO poo 



dz 



[z — wt, w) E (z + ws — wt, s) dw+ 



dgcc 
dw 



[z — wt, w) E {z + ws — wt, s) dw + R [z, t) (6.6) 



where we remark that the definition of R is: 



R{z,t) 



[fe [Voo [t; Z-Wt,w))- fe (w)] dw + 

+ I 5oo (^oo {t; z - wt, w) , Voo {t; z - wt, w)) dw- 
ds 

J —OO 

{z — wt, w) E {z + ws — wt, s) dw 



sds 



fe,v (w) E (z + WS — wt, s) dw — j goo {z — wt, w) dw+ 
' dgo 



dz 



OO /'OO 

ds 



dgoo 

dw 



[z — wt, w) E (z + ws — wt, s) dw 



(6.7) 



Using Assumptions (A), expanding the arguments of fe, goo using Taylor's 
expansion as well as the estimates (|6.4I) . (|6.5p we obtain: 



Riz,t) <C\\Ef{t + iye 



2 27* 



(6.8) 



We now transform the third term on the right side as suggested above, using 
integrations by parts. It then follows that (see p.l7p V 



Ez{z,t)= I ds fe.v (w) E (z — wt + ws, s) dw + / goo {z — wt,w) dw 

J t J — OO J —OO 

nOO ^ nOO 

+ / (t — s)-ds / goo{z — wt,w) Ez{z — wt + ws,s^ dw+ 

Jt t J-oo 

s\ dgo 



OO /"OO 



ds 



1 



1 1 dw 

A crucial estimate is the following: 



(2; — wt, w) E (z — wt + ws, s) dw + R (z, t) 



{t — s) -ds I goo {z — wt, w) Ez {z — wt + ws, s) dw 
t J —00 

poo 

< Ce \\E\\ e"^' / (s - t) -e'^^'-^^ds < Ce \\E\\ e^'^* 
Jt t 



(6.9) 
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where e comes from goo- A similar estimate can be obtained for the fourth term. 
The problem is then ready for a fixed point argument in the form 

/>oc> />oo 

Ez{z,t)— I ds / fe^jj (w) E {z — wt + ws,s) dw — {z,t)+L {z,t) + R{z,t) 

J t J — oo 

(6.10) 

where: 

/oo 
goo {z — wt,w) dw (6-11) 
-oo 

|.oo g |.oo 

L{z,t)— I {t — s) -ds I goo {z — wt,w) E2 {z — wt + ws, s) dw+ (6.12) 



t 



rOO ^00 ^ 

/ ds (1 I °° (z — wt, w) E (z ~ wt + ws, s) dw 

Jt J -00 ^ tl dw 

More precisely, the fixed point scheme is the following one. Given E with 
||i?|| < 00 we define L as in (|6.12p and R as in (|6.7p . We then define T {E) as 
the solution of (|6.10p obtained by means of Proposition [TBI 

In particular, using Lemma [SI (|6.8p . (|6.9p and Proposition [THj we have the 
following estimate: 

\\T{E)\\<Cg^e + Ce\\E\\+C\\E\\' (6.13) 

where Cg^ is a constant of order one that comes from g^o and depends only on 
B, a. 

A consequence of (|6.13l) is that the ball < 2Cg^e is transformed by 
means of the operator T in a set contained in a ball with radius: 

and this number is strictly less than 2C'g^e if e is sufficiently small. 

Finally we will prove that the operator T is contractive in the ball < 
2Cg^£ if e is sufficiently small. To this end we denote as Zi (s, t; z, v) , Vi (s, t; z, v) , 
Z2 (s, i; z, v) , V2 (s, t; z, v) the evolution of the characteristics defined by means 
of (ISH)- ([63]) with electric fields respectively with || <2Cg^e, \\E2\\ < 

2Cg^e. We will use also the notation Zi^^o {t; z, v) , Vi^oo {t; z, v) , ^2,00 {t; z, v) , 
"^2,00 (t; z, v) , Li (z, t) , L2 {z, t) , Ri (z, t) , R2 (z, t) to denote the corresponding 
functions associated with Ei, E2- 

Notice that ^^-^^ imply: 



d{Z^-Z2) 
ds 

d{Vi~V2) 
ds 



with: 



-s [El (Zi + V^s, s) - E2 {Z2 + V2S, s)] 
[El (Zi + Vis, s) - E2 {Z2 + V2S, s)] 



{Zi- Z2)it,t;z,v) = (6.14) 
{Vi-V2){t,t;z,v) = 
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We can write the equations as: 

d{Zi-Z2) 



ds 

d{Vi-V2) 
ds 

where: 



-s [a [s) (Zi -Z2)+b (s) {Vi - V2) + {El {Z2 + V2S, s) - E2 {Z2 + V2S, s))] 

(6.15) 

[a (s) (Zi -Z2)+h (s) (Fi - V2) + {El {Z2 + V2S, s) - E2 {Z2 + V2S, s))] 

(6.16) 



, , El (Zi + Vis, s) - El {Z2 + Vis, s) 

^^'^ ^ ^T^^ 

El {Z2 + Vis, s) - El {Z2 + V2S, s) 



b{s) 



V -V2 



Notice that: 



\a{s)\ + \b{s)\<\\Ei\\e-''^ (6.17) 
\Ei{Z2 + V2S,s)-E2{Z2 + V2S,s)\ < \\Ei~E2\\e-'" (6.18) 

Then, combining (|6.14p - (|6.18p we obtain, using a Gronwall type argument: 

\{Zi - Z2) {s, t- z,v)\ <C\\Ei- E2\\ {t + 1) e-T* 
\{Vi - V2) {s, t-z,v)\<C \\Ei - E2\\ e-T* 

whence: 

\{Zi,oo-Z2,oo){t;z,v)\<C\\Ei-E2\\{t+l)e-'^' , (6.19) 
\{Vi,oo - V2.,oo) {t; z, v)\ < C \\Ei - E2\\ e-^* 

Let us denote as Ei, E2 respectively the solutions of (|6.10p with sources 
Li, Ri and L2, R2 respectively. 
Notice that: 

(^Ei — £2^ {z, " y '^^ J ^'^'^ ^^-^ ~ i'^ ~ ^ "^^j ^) 

(6.20) 

= {Li~L2){z,t)+(Ri-R2) {z,t) 



Notice also that, arguing as in the estimates of L above (cf. (|6.12|) and 

my- 

\{Li-L2){z,t)\<Ce\\Ei-E2\\e-''' , zeM. , t>0 (6.21) 
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We also need to estimate the difference (-Ri - • Therefore we need to 
estimate: 

(Ri-R2]{z,t)^ / [fe {Vi^oo {t;z - wt,w)) ~ fe {¥2,00 {t; z - wt,w))]dw+ 

J — oo 



— 00 
00 



+ / [goo {Zi,oo {t; z - wt, w) , Vi,oo {t; z - wt, w)) - 
-goo (^2,00 (t; z - wt, w) , ¥2,00 [t; z - wt, w))] dw+ 

fe,v (w) {El [z + ws — wt, s) — E2 {z + ws — wt, s)) dw+ 
dgo 



00 roQ 
ds 



sds 



dz 



{z ~ wt, w) {El {z + ws — wt, s) — E2 {z + ws — wt, s)) dw- 



00 /'OO 



ds 



dgo 



dw 



{z — wt, w) {El {z + ws — wt, s) — E2 {z + ws — wt, s)) dw 

(6.22) 



Taking into account (|4.2|) and (|6.19p we can estimate the second, fourth and 
fifth terms in (jg?^ as Ce \\Ei - E2\\ e"^*. Then: 



(i?i-i?2) (z,t)| < J 



[fe {Vi,oo {t;z- wt,w)) - fe {¥2,00 {t; z - wt,w))] dW' 



00 pOO 

/ fe.v i'w) {El {z + ws — wt, s) — E2 {z + ws — wt, s)) dw 

'J —00 

+ Cs\\Ei- E2\\e-'^' (6.23) 
On the other hand (|6.2p . (|6.3|1 imply: 

/oo 
Ei {Z, (s) + sVi (s) ,s)ds , i^l,2 (6.24) 

where, for simplicity Zi (s) = Zi {s, t; z — wt, w) , Vi (s) = Vi {s, t;z — wt, w) , i = 
1,2. 

Using Taylor's Theorem and (|6.19p we can estimate the difference 
[fe (Vi,oo {t; z - wt, w)) - fe (V^2,oo {t', z - wt, w))] as: 

[fe {Vl^oo {t; Z - wt, w)) - fe (1^2,00 {t; Z - wt, w))] 

= fe,v {V2,oo {t; Z - wt, w)) {Vi^oo {t; Z - wt, w) - ¥2^00 {t] Z - wt,w)) + 

+ o(\\Ei~E2\\''e-^^'^ 



Using this estimate as well as (|6.19|) . (|6.24p and the fact that ||£'i|l < 
2Cg^e, i = 1, 2 we obtain: 

[fe {Vl^oo {t; Z - Wt,w)) - fe (^2,00 {t; Z - wt,w))] 

= fe,v {w) {Vi,oo {t; z - wt, w) - ¥2,00 {t; z - wt, w)) + 
+ O (\\Ei - i;2||' e-^T*) +0{e\\Ei- E2\\ e^^') 



32 



Taking into account (|6.24p it then follows that: 

[fe {Vi^oo {t;Z~ Wt,w)) - /e (V2,oo [t; Z - wt,w))] dw- 



< 



ds 



ds 



fe,v (w) {El {z + ws — wt, s) — E2{z + ws — wt, s)) dw 



OO 
OO 



fe,v (w) {[El {Zi (s) + sVi (s) ,s)^ Ei{z + ws- wt, s) 

) 

[E2 {Z2 {s) + sV2 {s) ,s)~ E2{z + ws- wt, s)])\ + 
C (\\Ei - ^2||' e-27t +e\\Ei- ^2!! e"^*^ 



Notice that 

{El {Zi (s) + sVi {s) , s) -Ei{z + ws~ wt, s)) 

- {E2 {Z2 (s) + SV2 (s) , s) -E2{z + ws- wt, s)) 

= {El {Zi {s) + sVi (s) ,s)-Ei {Z2 (s) + sV2 {s) , s)) + 
+ {El {Z2 (s) + SV2 (s) , s) - Ei{z + ws- wt, s)) 

- {E2 {Z2 (s) + sV2 (s) , s) -E2{z + ws~ wt, s)) 

The term {Ei {Zi {s) + sVi (s) ,s) — Ei {Z2 (s) + SV2 {s) , s)) can be estimated 



as: 



Zt{s)+sVi{s) 



El,. {Ls)di 



< C \\Ei\\ e-^' [\Zi (s) - Z2 (s)| + s 1^1 (s) - V2 (s) 



Z2{s) + sV2{s) 

whence, using \\Ei\\ < 2Cg^e: 

\Ei {Zi {s) + sVi {s),s) - El {Z2 (s) + sV2 {s),s)\ < Cee-^^e"-^* (1 + s) {\\Ei - E2 

(6.25) 

On the other hand 

{El {Z2 (s) + SV2 (s) ,s)~ Ei{z + ws- wt, s)) - 
- {E2 {Z2 (s) + sV2 (s) ,s)-E2{z + ws- wt, s)) 

^Z2{s)+sV2(s) 

[El^z {Ls)^E2^z {Ls)]d^ 



z+ws—wt 



thus we have 



I {El {Z2 (s) + SV2 (s) ,s)~Ei{z + ws~ wt, s)) - 
- {E2 {Z2 (s) + sV2 (s) ,s)~ E2{z + ws~ wt, s))\ 

< \\Ei ~ E2\\ e-^' [\Z2 (s) -{z + ws-wt)\+s \V2 {s) - w\] 

< Ce-^'e--" (1 + s) \\Ei - E2\\ \\E2\\ (6.26) 
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Then: 

[fe (Vl,oo (i; Z - Wt,w)) - fe (V2,oo [t] Z - wt, w))] dw- 



OO 

oo 



ds / fe,v i'w) {^1 {z + ws — wt, s) — E2 {z -\- ws — wt, s)) dw 

J — 00 

< Ce\\Ei- E2\\e-''* 
and therefore: 

(^Ri - R2) {z,t) <Cs\\Ei- E2\\e-^' (6.27) 
Combining dnHH), dOOl), fOTI and Proposition [H we obtain: 

\\T (El) - T {E2)\\ < Ce\\E,~E2\\ 

for \\Ei\\ < 2Cg^e, i — 1,2, thus the desired contractibihty of T foUows. This 
concludes the proof. ■ 
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